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We show how the noise in a spin polarized STM tunnehng current gives valuable spectroscopic 
information on the temporal susceptibility of a single magnetic atom residing on a non-magnetic 
surface. 
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I. INTRODUCTION 

No fundamental principle precludes the measurement 
of a single spin, and therefore the capability to make 
such a measurement depends on our ability to develop a 
detection method of sufficient spatial and temporal res- 
olution. The standard electron spin detection technique- 
electron spin resonance- is limited to a macroscopic num- 
ber (> 10^") of electron spins [1]. Recent experiments 
employing spin polarized STM [2] open new possibilities 
of investigation of magnetic systems at spatial resolutions 
of the Angstrom scale. The bulk of the experiments per- 
formed to date have been on magnetically ordered states, 
such as antifcrromagnets. 

Here we propose to use a spin polarized STM tun- 
nel current to gain spectroscopic information on a single 
magnetic atom on a non-magnetic surface. The set up is 
similar to the one used in recent ESR-STM experiments, 
although we now consider the same for spin polarized 
current. A scheme is depicted in Fig.(l). We start with 
a magnetic atom of spin S placed on an otherwise non- 
magnetic substrate. As wc will explain in the text, the 
noise in the current flowing from the STM tip will allow 
us, in certain instances, to directly measure the single 
spin time dependent susceptibility. This is yet another 
case where the spectroscopy of current noise allows us 
(where other methods often fail) to directly probe the 
highly disordered quantum states of a microscopic sys- 
tem (in this case a single spin). 

In earlier works, one of us and others examined simi- 
lar issues at finite magnetic fields. This work is crucially 
different from its predecessors in that we consider the sit- 
uation at a vanishing external field, and within the con- 
text of a spin polarized tunnel current. Our work also 
differs in the former context from [3], and further offers 
a transparant and accesible account of the underlying 
physics. Throughout the paper we propose and analyze 
aspects of a new experimental technique for probing sin- 
gle spin dynamics. 

The outline of this article is as follows: In section(II), 
we derive the general relation between the noise 



(|(5/(w)p) in the current and the single spin susceptibil- 
ity x{^)- Under quite general circumstances, the noise 
in the current is directly proportional to the single spin 
susceptibility, so that measuring the noise in the current 
immediately gives us valuable information about the sin- 
gle spin susceptibility. In Section(III), we will elaborate 
on the origin of the spin dependent tunneling matrix el- 
ements that form much of the backbone of our proposal. 
In Section(IV), we will discuss the decoherence result- 
ing from backaction effects. Here, we will also relate the 
spin scattering relaxation rate and DC charge transport 
and estimate how large the various quantities might be 
in realistic setups. We will conclude, in Section(V), by 
evaluating the typical signal to noise ratios for our ex- 
perimental proposal. 




FIG. 1. The experimental setup. The surface tip separa- 
tion distance is d. The single magnetic atom on the surface 
is marked by a solid circle. The spin polarized current ema- 
nating from the STM tip is schematically shown by the spin 
up open symbols. 
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II. HOW TO PROBE THE SPIN 
SUSCEPTIBILITY VIA A MEASUREMENT OF 
THE NOISE 

There is a major difference between probing a single 
spin (as discussed here) and probing a macroscopic mag- 
netic system. In the case of a single atomic spin S(i), 
we perform a measurement on a microscopic state that 
is quantum and highly fluctuating due to interactions 
with its environment. The tunneling experiment dis- 
cussed here is a specific example of noise spectroscopy, 
wherein we extract spectroscopic information, such as 
the single spin relaxation time T* from the measurement 
of the noise in the current. The main idea of noise spec- 
troscopy in STM is that the time dependence of the tun- 
neling current will have a characteristic relaxation time 
that is directly related to spin relaxation time T* . Here 
T* does not originate from a directly applied external 
magnetic field (there is none) but rather from substrate 
excitations and other effects. 

The time varying tunneling current I{t) will have a 
fluctuating component, apart from its average DC value 
Jq. Current noise is given by the Fourier transform of the 
time-dependent fluctuations of the electrical current lead- 
ing to the power spectra (|/(ti))p) within the frequency 
domain. Here, we will focus on the power spectrum of 
the tunneling current and argue that it is proportional 
to the local spin susceptibility, 

l/T* 

of the single atom on the substrate. In the above. T* de- 
notes the relaxation time to the polarization axis of the 
single local spin (for a macroscopic collection of spins, 
this would be none other than the standard Ti of NMR 
gauging the relaxation time to the polarization axis). We 
reiterate that in the present discussion there is a direct 
externally applied magnetic field T*, and albeit the sim- 
ilarity the origin of the relaxation time here and in NMR 
of a completely different origin. 

A long time exponential relaxation for x(t) = 
{S{t)S{0)) ~ exp[-t/T*], has as its Fourier transform 
the Lorentzian of Eqn.(l). In the Lorentzian form of 
Eqn.(l), the inverse relaxation time (T*)~^ trivially ap- 
pears as the linewidth in x(w). If additional short time 
spin modulations occur, these will augment x('^) by ad- 
ditional high frequency terms. 

In what follows, most of our focus will be on the long 
time (short frequency) behavior of the current noise. In 
thermal noise, the spectral density (|/(a>)p) is indepen- 
dent of frequency (e.g. being AksT / R at a temperature 
T for current flowing through a resistor of magnitude 
R). Shot noise, on the other hand, results from the dis- 
crete pulse-like character of electrical current. Its mag- 
nitude exhibits a white noise (frequency independent) 



spectrum that is proportional to the average DC current, 
(|/(a;)p) ~ ae/o (where the numerical constant a = 2 for 
a simple conductor), up to a certain cut-off frequency 
which is related to the time required for an electron to 
travel through the conductor. 

Examples of noise spectroscopy include the noise NQR 
measurements [5], noise in Faraday rotation [6] and, re- 
cently, noise spectroscopy of a local spin dynamics in 
STM [7]. The central feature present in all of these ex- 
amples is that the quantum system is not driven by ex- 
ternal flelds. Rather, the noise in the signal itself (e.g. 
thermal, shot, ...) sans any applied polarizing field allows 
us to extract spectroscopic information. Here we con- 
sider the excess noise produced by a single spin whose 
time-dependent quantum state we wish to probe. We 
suggest that over a certain parameter range, the excess 
noise generated by the single spin will be the dominant 
noise source in a single spin tunnel junction. 

To make matters concrete, consider the tunneling be- 
tween two contacts in the presence of a localized spin S. 
The Hamiltonian of this system assumes the form 

+ E ^LJ*0 + hS ■ &a0]CRn'0} + {L^ R). (2) 

nn'a/S 

In the above, &af3 is the Pauli matrix vector with ma- 
trix indices a and f3, the fermionic c\ncr,c[n<j ^^e 
annihilation and creation operators of electrons in the 
n-th eigenstate of the lead X = L,R with a = ±1 the 
(up/down) spin polarization label. The left lead (L) is 
the STM tip, and the right lead (R) refers to the surface. 
For infinite "free" leads, the eigenstates {|n)} simply cor- 
respond to the various plane wave states {|k)}. In real 
systems, all hopping matrices t will carry additional n, n' 
labels. However, to make the notation more compact we 
will dispense with these indices. The wavefunctions of 
our system are superpositions of the direct product states 

IVi) ® IV's) ® IV-ii) (3) 

- the direct product of the state of the left contact, the im- 
purity spin, and the right contact. The tunneling matrix 
i present in the second term of Eqn. (2) couples all of these 
different states. It has two contributions: the term pro- 
portional to to describes the spin independent tunneling 
while the term proportional to ti depicts the spin depen- 
dent contributions arising from the exchange interaction 
for electrons tunneling to the magnetic atom. Only the 
second term in Eqn. (2) will give rise to net current flow 
from the left to the right contact. In section(III), we will 
explain the origin and elaborate on the magnitude of the 
spin independent and dependent terms. A chemical po- 
tential shift (voltage drop) between the two bands {eni} 
and {cur} will lead to a DC current within the steady 
state. 
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Henceforth, we will assume that the tunneling elec- 
trons are partially spin polarized. There are several sit- 
uations where such interactions may materialize. In a 
ferromagnetically coated tip a chemical potential differ- 
ence 2((5/Ltcr) separates the two different spin polarizations: 
^Lna = ^Ln + f 5/1^ . Ferromagnetically ordered tips have 
proven to be very successful in the study of magnetic 
structures [2]. 

Very recently, an alternative approach using antiferro- 
magnetically coated tips with no ferromagnetic order has 
been used to produce spin polarized current [4]. These 
tips might have potential benefits as compared to ferro- 
magnetic tips. A ferromagnetic tip may produce a field 
of 0(1) Tesla at a separation of few angstroms from the 
surface. Ferromagnetic tips may therefore lead huge pre- 
cession frequencies that are difficult to mcasiire. In the 
case of an antiferromagnetic tip, there is a vanishing dipo- 
lar field. 

Both of these techniques may be used for the mea- 
surements proposed here. In what follows, we are not 
interested in a specific model for how the spin polarized 
current is generated. We define a parameter A that re- 
lates the spin polarized current to the net tunneling cur- 
rent. It is this parameter that will be determined by a 
particular microscopic model of the tip. Hereafter, we 
will treat A as a phenomenological parameter. 

As may be seen by examining the spin dependent con- 
tribution to t, the tunneling electrons exert torques on 
the localized spin S which lead to corrections to the spin 
dynamics [8]. To lowest order in ti, however, such effects 
are not present. Similarly, in what follows, we ignore 
the spin-flip interactions between the local moment and 
electrons in the substrate (i.e. we assume that the exper- 
imental temperature is higher than any relevant Kondo 
temperature (T > Tk))- Below the Kondo temperature, 
we may not ignore the spin-flip interactions of the lo- 
calized spin with the substrate electrons- the local spin 
susceptibility will be heavily influenced by these interac- 
tions. Here we will only consider T > Tk (a free impurity 
spin) for the sake of simplicity. For temperatures lower 
than the Kondo scale, the Kondo effect might manifest 
itself through interesting changes in the observed current 
noise that w(^ discuss here for the free spin case. 

To make our expressions slightly more concise, we 
will employ the Heisenberg representation and absorb 
the time dependence in all operators {O}, i.e. 0{t) = 
exp[ii7i]Os exp[— ii/f], with Os — 0{t = 0) the oper- 
ator in the Schroedinger representation which we now 
forsake. All finite temperature expectation values {0{t)) 
that appear will represent X]iPi(V'i(0)|O(t)|V'i(0)) with 
Tpiit = 0) the zero time wavefunction of the Schroedinger 
representation which does not evolve within the Heisen- 
berg formulation and Pi its probability within the density 
matrix formulation. 

Let us now give a simple qualitative description of 
the effect that we address here. By directly computing 



dN^/dt we find that the charge current [9] 

I{t) = -ie ^ 4na W[*0^a/3 + tiS{t) ■ aa0]cRn'0{t) 

nn'ajS 

+h.c., (4) 

with e the electronic charge. 

Thus, we see that the tunneling current has a part that 
depends on the localized spin via a scalar product, 

6I{t)=etiS{t)-ispin{t), (5) 

where 

hpin{t) = -i ^ c{^ai*)^c,0CRn'f3{t) + h.C, (6) 
nn'af3 

is the spin polarization dependent contribution to the 
electronic current. This expression has a very transpar- 
ent meaning. Its z-component, /fpi„(i) = ^'-(c^ ^cr.'\ — 

^cr^i) + h.c, is the net flow of up spin minus the flow 
of down spin. A DC current of polarized electrons flow- 
ing from the tip to the surface trivially leads to a uni- 
form spin polarized current Ispin- We assume that there 
is a non-vanishing steady spin polarized current compo- 
nent tunneling from the tip to the surface, assumed to be 
aligned along (or deflning) the z axis: {Ilpin{i)) = Si,z^^ 
+ time dependent fluctuations, with a finite A ^ 0. An 
application of the diagrammatic analysis introduced in 
[10], reveals that four contributions result each of which 
is, at most, of the order of the contribution that we dis- 
cuss below. To lowest order in the hopping amplitude ti , 
the electronic current-current correlation function orig- 
inating from the spin dependent part that we wish to 
probe, 

{{sm,si{t')}) = eHl{s\t)s^{t')){il^^{t)il,^^^^^^ 

+{t^t'), (7) 

where {} denotes a symmetrized correlator, i,j = x,y,z 
denote the spin components. To lowest non-trivial or- 
der in ti, we need to treat the two temporal correla- 
tion functions xit - t'\ = (5^(t)S'^(t')) and C{t - t') = 

{ilpur(t)ilpM ^ {ilp^n{t)){iip^n{t')) = k^h^A (for 

\t — t'\ oo) independently. To this order, the wave- 
functions with respect to which we compute the ex- 
pectation values are the products of pieces describing 
the decoupled evolution of both the spin and of the 
left and right contacts separately. To make connection 
with the main proposal in this paper, we note that, 
when Fourier transformed, the symmetrized correlator 
{{5I{t),5I{t'y}) is none other than the current noise spec- 
trum at various frequencies originating from the local 
spin. For small (ii/to) ^ 1 (the experimental situation), 
0{{{51{t),5i{t')})) = {ti/tofll. In evaluating the spin 
current correlator C{t), we ignore the fluctuating contri- 
butions present for short times (large frequencies). The 
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spin current correlator C{t) has a finite, asymptotic, long 
time value. A, that reflects the spin polarized DC current 
emanating from the STM tip. In Fourier space, the cur- 
rent power spectrum is given by a convolution of the two 
power spectra associated with S (i.e. x{^)) a.nd cr (the 
spin current correlator C{co)): 



2 2 f dUi 



LiJi)C(0J — (jJi) + (OJ 



(8) 



At low frequencies, C{w) ~ 2'kA5{w), and, consequently, 
{\5i{u^)\'')=2AeH\x{<^) + ... (9) 

The ellipsis in Eqn.(9) refer to the contribution to the 

convolution of Eqn.(8) from the finite frequency (short 
time) contributions to C{t). Assuming such finite fre- 
quency contributions in C{uj) have low spectral weight, 
the effect of these contributions will be low. This low or- 
der result is augmented by higher order corrections in ti 
as well as contributions arising from the connected com- 
ponent of the current current correlator which amounts 
to a shot noise like contribution. 

As noted, in Eqn.(9) we neglected the effect of fi- 
nite frequency components of C{uj) in the convolution 
of Eqn.(8). We beheve the large finite frequency com- 
ponents of the spin current C{w) to be small as these 
correspond to transient fluctuations about an assumed 
steady state. This assumption is not necessary, however. 
More generally, we may deal with the full convolution in 
Eqn.(8) directly without invoking any assumptions. We 
may potentially do this by first performing a measure- 
ment on a reference state. For instance, if we initially 
replace the single spin by a large cluster of a large fixed 
spin, we may then experimentally measure the resulting 
spin current C{w). Armed with the knowledge of the spin 
current (from this earlier measurement on the magnetic 
cluster), once the noise spectrum (|/(a;)p) is measured 
for the single spin, we can directly deconvolve Eqn.(8) to 
obtain the spin susceptibility [11]. 

Eqs.(8, 9) are our central result. They vividly illustrate 
how the spectroscopy of the noise in the tunneling cur- 
rent (|(5/(ti))p) allows us to directly probe the spectrum 
of spin fluctuations encapsulated in x(a;). The spin po- 
larized tunneling current provides a reference frame with 
respect to which wc may measure the fluctuations of the 
localized spin S(i). We note that these results are gener- 
ally applicable to systems displaying strong correlations 
(such as impurities in the Kondo regime). 



III. THE ORIGIN OF SPIN DEPENDENT 
TUNNELING 

We now elaborate on the origin and magnitude of the 
spin dependent tunneling matrix elements of Eqn.(2). 



The spin dependence of the tunneling originates from the 
direct exchange dependence of the tunneling barrier [7] . 
The overlap of the electronic wave functions of the tip 
and surface, separated by a distance d is exponentially 
small and is given by a spin dependent tunneling matrix 
element, 



i = 7exp[ 



JS(t)-<T, 



$0 



(10) 



where we explicitly include the direct exchange between 
tunneling electron spin cr and the local spin S. Here, J is 
the exchange interaction between the electrons tunneling 
from the tip and the local precessing spin S. In the above, 
i is to be understood as a matrix in the internal spin 
indices, and <I> is the tunneling barrier height. Typically, 
$ is a few eV. As a canonical value we may assume $ = 
4eV, $0 = smP related to the distance d between 
the tip and the surface [12]. As the exchange term in 
the exponent is small compared to the barrier height, we 
may expand the exponent in JS. Explicitly, t may be 
written as 



i=to + ti&-S{t), 



where. 



to = 7exp(-($/$o)^^^)cosh[ 




(11) 



(12) 



describes spin independent tunneling. The spin depen- 
dent amplitude. 



ti =7exp(-($/$o)^/^)sinh[ 



(13) 



'2$ V $0 

For estimates we may employ the typical rule of thumb 

JS 
2* 



h/to ^ il « 1- 



IV. BACKACTION EFFECT OF THE 
TUNNELING CURRENT ON THE SPIN 

We may use the timncling Hamiltonian of Eqn.(2) to 
estimate the decay rate of the localized spin state result- 
ing from the spin scattering interaction associated with 
ti. To second order this calculation is equivalent to a sim- 
ple application Fermi's golden rule leading to an up-down 
spin flip rate = TrtlNLNjieV, with V is the voltage 
applied between the left and right electrodes. Similarly, 
the DC tunneling current /q is given by the tunneling 
rate of conduction electrons ^ = nt^Ni^NfieV, where 
^L,R denotes the density of states at the Fermi level of 
the tip and surface respectively [13]. 

Diagrammatically, both the spin scattering and elec- 
tronic scattering rates arise from the simple bubble 
diagram whose real-time propagators are GL{t) and 
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Gii{—t), where Gl.b. correspond to the left and right 
Green's functions for the conduction electrons respec- 
tively. The sole difference between the two (spin depen- 
dent/independent) scattering rates is encapsulated in the 
prefactors. In the spin dependent scattering case (t~^) 
the raw value of the single loop integral appearing in the 
bubble diagram needs to be scaled by tf. The spin in- 
dependent scattering rate (r~^) is the much same albeit 
a scaling by the spin independent scattering amplitude 
squared (^q). 

Comparing these simple results, we find 



Ts to Te 2$ 



(14) 



The important outcome of this analysis is that the cur- 
rent induced broadening predicts a spin relaxation rate 
— cx /q which may be experimentally tested. This re- 
sult has a very simple interpretation: the impinging for- 
eign electron tunneling rate for a DC current of magni- 
tude Iq = InA is given by 10^"i/z. By contrast, 
the probability to produce a spin flip, sparked by the 
tunneling electrons, is proportional to i^, which leads to 
Eqn.(14) for the linewidth. The full intrinsic line width 
is further enhanced by the coupling of the spin to the 
environment (e.g. the interaction between the spin and 
various substrate excitations) which may indeed further 
increase the spin flip rate. The net observed linewidth, 

{T*)-'^r-'+r-„l, (15) 

includes both backaction contributions (t~^) and the 
aforementioned linewidth broadening due to coupling to 
the environment (r^^^). The inverse backaction relax- 
ation time scale sets a lower bound on the net relaxational 
linewidth of the single impurity spin. 

Given the typical values of the parameters in Eqn.(14), 
we estimate ~ 5 x lO^Hz for our DC current of 
lo InA, J - 1 eV, $ - 4 eV, and S = 1/2. Future 
experiments will help to clarify the linewidth dependence 
on the various parameters. 



V. A SIZABLE SIGNAL TO NOISE RATIO 

We now demonstrate that the (signal to noise) ratio 
of \5{oj — > 0)p (the excess noise induced by the impurity 
spin) to \Ishot{i^ ^ 0)p (the shot noise already present 
in the absence of the impurity spin) can be of quite sig- 
nificant (of order unity). 

The finite frequency ratio \6I{co)\'^ /\Ishot{'^)\'^ was 
computed in a multitude of systems and shown to be 
bounded by 4 or other numbers of order unity. By con- 
trast, the low or zero frequency signal to noise ratio 
\61{u! O)p/|/s?iot(w 0)|2 was found to be unbounded 
in many instances. For one calculation amogst many oth- 
ers demonstrating this explicitly, see e.g. [14]. Neverthe- 
less, as we highlight below, the low frequency signal to 



noise ratio, \5I{lj O)]"^ /\Ishot{^ 0)1^, in our sys- 
tem is stringently bounded from above by a number of 
order unity, just as it is in many finite frequency situa- 
tions. By considering typical parameter values, we will 
show that this ratio may, potentially, saturate this up- 
per bound and be markedly large. The spin dynamics as 
manifested through the current noise may be very easily 
discernible. 

To estimate this ratio, we note that the average spin 
dependent contribution to the current. 



{6m=et,{S\t)P^Jt)) 



(16) 



leading to 0{etiA^/'^{S''{t))). Here, as throughout, H = 
1, and the spin is dimcnsionlcss. 

To obtain estimates of orders of magnitude let us in- 
spect Eqs.(4, 6). The net electronic current Ig of Eqn.(4) 
is of the order of Iq. Insofar as orders of magnitude arc 
concerned, the spin current {I spin) defined in Eqn.(6) sat- 
isfies 



The shot noise, (Ighot) = ^eJo, where the numerical 
constant a = 0(1). Noting that Iq = and making use 
of Eqn.(9), the signal to noise ratio is found to be 



(|«(a;)n 



In 



1 

Te 



(18) 



The ratio {h/to) — 55. Inserting the Lorentzian form of 
X(a;) from Eqn.(l), we obtain 



(|5W(a;^0)|2) 



T* ,t 



1 \2 



if) 

J rp^ 
{ — ) — 



(19) 



Inserting typical values, J ~ 0.1 — 1 eV, $ ~ 4 eV, 

T* - 10"*-10~^ seconds, andr,. - 10-^° seconds (InA), 
we find the above signal to noise ratio is, naively, 1-100. 
As promised, we now demonstrate, within this allowed 
empirical regime, the signal to noise ratio will typically 
veer towards the lower end of the spectrum (i.e. may be 
of order unity at most). The bottleneck in the signal to 
noise ratio is set by the backaction effect of the single 
spin on the tunneling current. More explicitly, fusing 
Eqs.(14,15) together, 



(J^*)-l ^ 

\ / — ' s ' ' env 

~ (Z:^)2l , -1 

leading us to conclude that T* < {^)'^Te. 

Inserting this back in the last line of Eqn.(19), 

{\6hhot{i^^0)\^) -4$^ Te ~ ^ 



(20) 



(21) 
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As the generic signal to noise ratio is bounded both 
from below and above by a number of order unity, the 
signal may indeed be quite sizable (of the order of the 
shot noise) and may be detected. 

This sizable ratio offers promise for such single spin 
detection and related small system applications. 

Note added. While finishing this paper we became 
aware of work by Bulaevskii, Hrushka and Ortiz [3], 
where a similar problem was considered for a specific case 
of a ferromagnetic STM tip. Our results are qualitatively 
similar to the results obtained by Bulaevskii et al. 
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